Transverse localization in nonlinear photonic lattices with second-order coupling 
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We investigate numerically the effect of long-range interaction on the transverse localization of 
light. To this end, nonlinear zigzag optical waveguide lattices are applied, which allows precise tuning 
of the second-order coupling. We find that localization is hindered by coupling between next-nearest 
lattice sites. Additionally, (focusing) nonlinearity facilitates localization with increasing disorder, as 
long as the nonlinearity is sufficiently weak. However, for strong nonlinearities, increasing disorder 
results in weaker localization. The threshold nonlinearity, above which this anomalous result is 
observed grows with increasing second-order coupling. 



I. INTRODUCTION 

The localization of waves in disordered lattice systems, 
as proposed in 1958 by P.W. Anderson [1], is a fascinating 
feature of energy transport. It was predicted that, when 
disorder is introduced to a periodic system, the extended 
Bloch eigenmodes may convert to exponentially localized 
states, leading to metal-insulator transition [2|. The in- 
terference between multiple scattered electronic waves is 
the origin of Anderson localization. Although, the An- 
derson localization was initially introduced in solid-state 
physics, due to its wave nature based on interference ef- 
fects, this concept could be applied to other waves such 
as light [sl-lil , matter [6] , and even sound . Notably, all 
of these studies do not consider the impact of non- nearest 
lattices sites on localization. However, in some systems 
like biomolecules [8] and polymer chains [9], long-range 
interactions in the lattice become important and cannot 
be neglected. 

In recent years, coupled optical waveguide lattices pro- 
vide an excellent platform for study of transverse local- 
ization of light [El, [lo|, [Hi. In particular, optical lattices 
fabricated using the femtosecond (FS) laser writing tech- 
nology are particularly useful to study various effects as- 
sociated with disorder [12-16]. Besides these experimen- 
tal studies of light localization in disordered photonic lat- 
tices, there is a large number of numerical or theoretical 
literature on the different aspects of these systems pjf . 
For instance, the competition between disorder and non- 
linearity [18|, dependence of localization on input beam 
profile p^, t ransverse localization with dimensionality 
crossover and the effect of the excited site number 
on the wave-packet localization have been studied 
previously. 

In our work, we study numerically the impact of in- 
teractions between non-nearest lattice sites on Anderson 
localization in a disordered system in the linear and non- 
linear regime. To this end we employ zigzag lattices of 
optical waveguides [22-25], which are particularly useful 
to study second-order coupling (SOC) between the lattice 
sites. 



This paper is organized in four sections. Section II con- 
tains the theoretical considerations of disordered zigzag 
arrays of optical waveguides, and the definition of stan- 
dard quantities for investigation of transverse localization 
in these systems. In section HI, the results of numerical 
simulations and discussions will be presented, and finally, 
conclusions will be presented in section IV. 



II. THEORETICAL CONSIDERATIONS 

Our system consists of a zigzag array of identical single- 
mode circular optical waveguides (Fig. 1). We assume 
that the angles between adjacent arms are the same (^), 
while the distance between two successive waveguides 
are selected randomly from an uncorrelated exponential 
probability distribution. According to the exponential 
dependence of coupling coefficient with waveguide sepa- 
ration [26], this will lead to a uniform random distribu- 
tion of coupling constants between adjacent guides. 

In the coupled mode approximation, the evolution 
equation for the complex amplitude of electric field at 
the jth waveguide Ej is |22|] 



+t^{j)Ej_2 + x\Ej\^Ej 



(1) 



with j = 1, . . . , A/", where tf{j) are the first-order cou- 
pling (FOG) coefficients between guides j and j±l, tf{j) 
are the SOC coefficients between non- nearest neighbors 
(NNN) j and j ± 2, and x is the nonlinear Kerr constant 
(we consider only focusing nonlinearities x ^ 0)- 

In lossless systems t+(j) = t~(j + m),m = 1,2 [2^. 
The FOC coefficients tf{j) are random numbers with 
uniform distribution in the interval [to{l — A), to(l + A)], 
where < A < 1 is the disorder strength. We set the av- 
erage distance between successive waveguides {duj) to 18 
/im, corresponding to to 0.15 mm~^ (for A = 780 nm) 
[2^. For each realization {dj}^ the SOC coefficients t^(j) 
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FIG. 1. Zigzag array of optical waveguides. The distance 
between guides j and j+1 is denoted by dj, and the angle 
between adjacent arms is 0. 
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can be determined by the NNN separation [26 
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2djdj-^i cosO 



(2) 



The relative coupling strength a, which is defined as 
/tf at A = 0, can be controlled with the angle in the 
Fig. 1. Not that at ^ < 60° one finds a > 1. The nonlin- 
earity in the system can be tuned by the Kerr parameter 
X, and X = corresponds to purely linear dynamics. 

To investigate the effect of nonlinear ity and SOC on 
Anderson localization, we solve Eqs . ([p numerically by 
the Runge-Kutta-Fehlberg method [23, with single site 
excitation {Ej{z = 0) = Sjj^) as initial condition. To 
consider finite size system effects, we use fixed boundary 
conditions, i.e., Ej = for j = —1, 0, + 1 and + 2. 
Under this condition the norm ^f^i \Ej\'^ is conserved. 
To peruse transverse localization, we use two measures: 
the dimensionless transverse localization length (TLL) 
which is obtained by fitting an exponential function 
\Ej\'^ = |E|^^^e~l-^~-^^l/^ on the localized profile, and the 
the effective width Weff [5]: 
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2^7 = 1 



that illustrates the formation of a localized state during 
the propagation. Due to the statistical nature of Ander- 
son localization, both quantities are averaged over 5000 
different realizations. 



III. RESULTS AND DISCUSSIONS 
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FIG. 2. (color online) Inverse of transverse localization length 
versus relative coupling strength a for different values of 
the nonlinear parameter x (mm~^). The excitation is at site 
jo = 50. (a) A = 0.4, (b) A = 0.7. 
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FIG. 3. (color online) Autocorrelation function of FOC and 
SOC coefficients for A - 0.7. 



For our simulation, we consider zigzag arrays contain- 
ing A" = 100 waveguides each has 100 mm length. In 
Fig. 2, we show the inverse TLL at the end of 

the waveguide array as a function of a for two disorder 
strengths A = 0.4,0.7 and several nonlinear parameters, 
when the lattice bulk was excited (jo = 50). The main 
result of our simulations is that, in general, SOC impedes 
localization and increases the transverse spreading of the 
wave function. According to the band random matrix 
theory [28|, it is well known that the localization length 
of the eigen states, at large values of the band width, in- 
creases with the square of the band width. The increase 



of the transverse spreading can be explained in terms 
of a modified band structure of the Hamiltonian matrix, 
in the presence of SOC. In the presence of coupling to 
the second- nearest neighbor, the Hamiltonian matrix of 
the system will be changed from tridiagonal form to a 
pentadiagonal one, hence the direct coupling paths for 
spreading of the wave packet will be increased. More- 
over, short-range correlations between the lattice sites, 
introduced by SOC, is another important feature of this 
system which may affect on the system localization [29| . 
In Fig. 3, we plot the autocorrelation functions of the 
FOC and SOC at A = 0.7, which indicates that although 
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FIG. 4. (color online) Inverse of transverse localization length 
versus relative coupling strength a for two case of corre- 
lated and uncorrelated SOC coefficients in the weak disorder 
regime A = 0.4. Excited site number is jo = 50, and (a) 
X = 0.1 mm~^ , (b) x — ^ mm~^ . 



FIG. 5. (color online) Inverse of transverse localization length 
versus relative coupling strength a for two case of corre- 
lated and uncorrelated SOC coefficients in the strong disor- 
der regime A = 0.9. Excited site number is jo — 50, and (a) 
X = 0.1 mm~^, (b) X = 1 mm~^ . 



the FOG coefficients are chosen from a uniform uncorre- 
lated random sequence, the SOC coefficients have posi- 
tive short-range correlation. This is due to the presence 
of the joint arm dj-^i in Eq. [2] for the expression of t2{j) 
and t2"(j + l). We find that the autocorrelation functions 
do not show any noticeable change with the variation of 
a. However, by decreasing the disorder strength A the 
value at Aj = 1 increases slightly. 

To investigate the effect of correlation on the localiza- 
tion length, we repeat our simulations but with uncorre- 
lated random hopping terms. It is important to notice 
that these uncorrelated random terms are chosen from 
the same probability distributions related to the model 
under study. Figs. 4 and 5 show the inverse TLL 
at the end of waveguide array as a function of a for two 
cases of correlated and uncorrelated SOC coefficients, for 
two disorder strengths A = 0.4, 0.9 and nonlinear pa- 
rameters X = 0-1 mm~^ and X = 1 mm~^. We find 
that site correlation effects appear only at high disorder 
strengths for small nonlinear parameters. In addition, in 
these regimes, correlation will decrease TLL {£) slightly 

Two selected propagation images for small {a = 0.3) 
and high {a = 2) SOC strengths, shown in Fig. 6, illus- 
trate our findings. The plots exphcitly demonstrate the 
enhanced expansion of the propagating wave packet for 
higher a, i.e., with increasing SOC. 

Another important result of our simulations is the en- 
hancement of localization by increasing nonlinearity (see 
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FIG. 6. (color online) Propagation images for A = 0.7, x — 
1 mm~^, jo = 50 and (a) a = 0.3, (b) a = 2. 



Fig. 7). In fact, one identifies two distinct nonlinear 
regimes, where the inverse TLL has small and large 
values [21], respectively. When the nonlinear parame- 
ter X is below some characteristic value Xc (see Fig. 7) 
(weakly nonlinear regime), localization growth slightly 
with increasing nonlinearity. In contrast, when the non- 
linear parameters is larger than the critical value X > Xc 
(strongly nonlinear regime), self-trapping effect arises 
and results in the formation of highly localized modes 
[U, [32|. Moreover, Fig. 7 shows that this characteristic 
value Xc increases with increasing SOC strength. 

It is interesting to analyze the impact of the disorder 
strength A on the transverse localization in the case of 
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FIG. 7. (color online) Inverse of transverse localization length 
versus nonlinear parameter x two values of relative 
coupling strength a. The excitation is at site jo = 50, and 
A = 0.7. 



non- vanishing nonlinearity. This topic already attracted 
much attention in recent years [33|, [s^ • In the absence of 
SOC, for X ^ 0-6 mm~^ (which is approximately equal 
to Xc] see Fig. 7 for a = 0) localization is stronger for in- 
creasing disorder level. This result is consistent with our 
common sense, because the interference due to the dis- 
order is an essential resource for Anderson localization. 
However, for x > Xc an increasing disorder level leads to 
less localization of the propagating wave packet. We at- 
tribute this anomalous result to destructive interference 
effects due to the nonlinear Kerr phase shift. In the pres- 
ence of SOC, for small values of a where localization is 
not hindered, a similar behavior is observed. However, in 
this latter case, the anomalous result appears for larger 
value of X, which is due to the growth of Xc with increas- 
ing relative coupling strength a. These findings are sum- 
marized in Fig. 8, where the effective width Weff of the 
wave packet is plotted as a function of at relative SOC 
strengths a = and a = 0.5, nonlinearity x = 0.7 mm~^ ^ 
and different disorder levels A = 0.4,0.5,0.6,0.7. One 
can clearly see that for such high nonlinear it ies, the pres- 
ence of SOC causes a decrease of the effective width with 
increasing disorder level. If we increase in our simula- 
tions the nonlinearity to x = 1 mm~^ - which is larger 
than Xc for both a = and a = 0.5 - the anomalous 
result appears again. In this case, for both a = and 
a = 0.5, high disorder levels lead to less localized wave 
packets. 

In order to complete our analysis, we repeat our anal- 
ysis on the inverse TLL as a function a and x foi" an 
excitation at the edge of the lattice (jo = !)• The re- 
sults are summarized in Fig. 9. We find all features of 
the bulk excitation (compare Fig. 2) also for the edge 
excitation. The only difference between the cases of edge 
and bulk excitation is that the inverse TLL for the edge 
excitation is a highly non-monotonic oscillating function 
of the relative SOC strength a, for each nonlinear Kerr 
parameter x (see Fig. 9). Though, the general trends 
prevails that localization gets weaker as a increases. 
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FIG. 8. (color online) The (dimensionless) effective width of 
wave packet as a function of z for different disorder levels, 
with parameters jo = 50, x = O-'^ mm~^ and (a) a = 0, (b) 
a 0.5. 
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FIG. 9. (color online) Inverse of transverse localization length 
versus relative coupling strength a for different values of 
the nonlinear parameter x {mm~^). The excitation is at site 
jo = 1. (a) A = 0.4, (b) A = 0.7. 



In the linear regime (x = 0) and in the absence of SOC 
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(a = 0), the ratio of for edge excitation (jo = 1) to 
the value for bulk excitation (jo = 50) is approximately 
equal to 0.6 and 0.9, for disorder strengths A = 0.4 and 
A = 0.7, respectively. This confirms that, in the linear 
regime, surface modes are less localized than bulk modes 
due to boundary repulsion [l2|, [35| . This ratio increases 
by 50% when the disorder level A grows from 0.4 to 0.7. 
Hence, for strong disorder levels, the localization of sur- 
face and bulk modes is comparable [2l|. However, as can 
be seen in Figs. 2 and 9, nonlinearity can reverse this 
effect. In this case, bulk modes become more extended 
than surface modes, which is in agreement with recent 
results 

IV. SUMMARY AND CONCLUSION 

We have studied the effect of SOC on transverse local- 
ization in zigzag arrays of circular optical waveguides in 
the presence of Kerr nonlinearity. Our simulations reveal 
that increasing next-nearest neighbor interaction hinders 
localization. We attribute this behavior to a modified 



band structure of the Hamiltonian matrix in the pres- 
ence of the second-order coupling. For bulk excitation 
the dependence of inverse TLL is a monotonic function 
of the SOC, where for edge excitation the dependency is 
highly non-monotonic and oscillates for increasing SOC. 
Moreover, in the absence of higher-order interactions, in 
the strongly nonlinear regime (i.e., when the nonlinear 
parameter is larger than some critical value), counter- 
intuitively localization will decrease with increasing dis- 
order level. The critical value, above which this anoma- 
lous behavior happens increases with growing SOC. 
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